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The low temperature electrical conductance through correlated quantum dots provides a sensitive 
probe of the physics (e.g., of Fermi-liquid versus non-Fermi-liquid behavior) of such systems. Here, 
we investigate the role of level asymmetry (gate voltage) and local Coulomb repulsion (charging en- 
ergy) on the low-temperature and low-field scaling properties of the linear conductance of a quantum 
dot described by the single-level Anderson impurity model. We use the numerical renormalization 
group to quantify the regime of gate voltages and charging energies where universal Kondo scal- 
ing may be observed and also quantify the deviations from this universal behavior with increasing 
gate voltage away from the Kondo regime and with decreasing charging energy. We also compare 
our results with those from a recently developed method for linear and non-linear transport, which 
is based on renormalized perturbation theory using dual fermions, finding excellent agreement at 
particle-hole symmetry and for all charging energies and reasonable agreement at small finite level 
asymmetry. Our results could be a useful guide for detailed experiments on conductance scaling in 
semiconductor and molecular quantum dots exhibiting the Kondo effect. 

PACS numbers: 75.20.Hr, 71.27.H-a, 72.15.Qm, 73.63.Kv 
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I. INTRODUCTION 

Artificial nanostructures, such as semiconductor quan- 
tum dots, 1-1 magnetic atoms adsorbed on surfaces, 5 ~ 
and, molecules attached to leads, s ~ provide new realiza- 
tions of the Kondo effect of a local spin interacting anti- 
ferromagnetically with conduction electrons. In contrast 
to their bulk counterparts, these systems are also highly 
tunable, for example, via application of gate voltages to 
modify the energy levels of the quantum dot or molecule, 
or to tune the tunnel couplings between the leads and 
the dot. In addition, application of a finite transport 
voltage allows an experimental investigation of the ef- 
fects of strong correlations on non-equilibrium transport 
through these model nanosystems, thereby motivating 
also the development of new theoretical approaches for 
non-equilibrium . l 

Motivated by recent experiments on conductance scal- 
ing in correlated quantum dots exhibiting the Kondo 
effect, 1 ' 1 "'"'' we present in this paper a detailed study 
of the low-temperature and low-field scaling properties 
of the linear conductance of a quantum dot described 
by the single-level Anderson impurity model. Scaling in 
physical properties is a hallmark of the Kondo effect.' ' 
Thus, a Kondo model description of a quantum dot im- 
plies that the conductance G(T, B) is a universal function 
of T/Tq and g(i B B/k B To over all temperatures T and 
magnetic fields B, with microscopic parameters (such as 
the Kondo exchange J) only entering through the dy- 
namically generated low energy scale Tq (to be defined 
explicitly in Sec. Ill), with g, /ib, k B denoting the g- 
factor, Bohr magneton and Boltzmann's constant respec- 
tively. In particular for T < To or gUftB <C k B To the 



functions G(T,B = 0) = G(0,0)(1 - c T (T/T ) 2 ) and 
G(T = 0,B) = G(0,0)(1 - c B (g/i B B/k B T ) 2 ) exhibit 
Fermi liquid corrections about the unitary conductance 
G(0,0) with deviations which are universal in the sense 
that the coefficients ct — 7r 4 /16 and cb — 7r 2 /16 are inde- 
pendent of microscopic details. Actual quantum dot de- 
vices, however, have a finite charging energy, and they are 
more realistically described by an Anderson model. The 
finite charging energy, and the ability to change the level 
energy of the quantum dot with a gate voltage, allow for 
charge fluctuations (even in the "Kondo regime" of the 
quantum dot) and can give rise to deviations from the 
expected Kondo scaling. It is therefore of some interest 
to quantify the effect of increasing charge fluctuations on 
the values of ct and cb- Recently, this issue has also been 
addressed in Ref. 25 by using a renormalized perturba- 
tion theory on the Keldysh contour 16,28 formulated using 
dual fermions. 2 ' 1 . This approach, denoted henceforth 
as superperturbation theory (SPT), yields both the linear 
and non-linear conductance. In this paper we shall com- 
pare the predictions of this theory for the linear conduc- 
tance with results obtained within the numerical renor- 
malization group (NRG) approach. 32 

The outline of the paper is as follows. Section II 
describes the quantum dot model. Section III gives a 
brief description of the calculation of the finite tempera- 
ture linear conductance G(T,B) of the Anderson model 
within the NRG following the procedure in Ref. 35 . In 
Sec. IV, some Fermi liquid results for ct and cb are given, 
and in Sec. V we outline the SPT calculations of ct and 
cb, with which we shall compare. In Sec. VI we present 
results for the dependence of the coefficients cxand c B 
on charging energy and gate voltage (local level energy) . 




FIG. 1. (Color online) A strongly correlated quantum dot 
with charging energy U 2> A and level energy Ed connected 
to leads via tunnel barriers. The gate voltage V g ~ Ed allows 
changing occupation of the dot nd from nd = 1 for Ed = — U/2 
to nd = through a mixed valence regime with nd ~ 0.5 for 
Ed ~ 0. In the singly occupied configuration, shown here 
for Ed ~ —U/2, the dot has a well defined spin 1/2 and the 
Coulomb blockade excitations at Ed and Ed + U correspond to 
removing or adding an electron. The coupling of the spin 1/2 
to the leads results in the Kondo effect, which is manifested by 
the appearance of an additional many-body Kondo resonance 
at the Fermi level ef = at low temperatures T < To . This 
resonance is also reflected as a zero bias anomaly in the non- 
linear conductance dl/dV in experiments. 58 



The latter are compared with the corresponding results 
from SPT. We conclude in Sec. VII with a discussion of 
the relevance of our results for experiments on quantum 
dots. In Appendix A we give an alternative derivation 
of the discretization scheme of Campo et al. in Ref. 36, 
which we have used in the NRG calculations reported in 
this paper. This derivation is carried out for an energy 
dependent hybridization function following the procedure 
in Refs. and . In Appendix B, we provide details of 
the SPT calculation of cb in terms of renormalized pa- 
rameters. The relation between bare and renormalized 
parameters, required for comparing SPT results for cb 
and ct with the NRG results, is also described in Ap- 
pendix B. 
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FIG. 2. (Color online) Linear conductance G(T) versus 
T/Tg for U/A = 16 and several values of s d = -U/2, 0, +(7/2 
using the approach of Ref. 35, with Tq defined by Eqs.(6-7). 
We also indicate with horizonatl and vertical dashed lines the 
extraction of the experimental Kondo scale, T£ xpt , from the 
mid-valley Kondo conductance (i.e. that for Ed — —U/2) via 
G(T = T™ pt ) = G(0)/2. NRG parameters were for A = 4, 
with an energy cut-off e c (A = 4) = 30 and n z = 2. 



Here, Ed is the level energy, related to the gate volt- 
age V g in the quantum dot via Ed ~ eV g (see Fig. 1), 
B is a local magnetic field acting on the quantum dot 
with sf = \{nd\ — nd\), U > is the Coulomb charg- 
ing energy, a labels the spin, and a = L, R labels left 
and right electron lead states with kinetic energies €fc Q . 
The couplings of the dot to the leads are denoted by 
A a (w) = irp a (uj)\t a \ 2 , where p a (oj) = J^k^i^ ~ e ka) is 
the density of states of lead a. For simplicity we as- 
sume a constant density of states p a = Np = 1/2D with 
half-bandwidth D = 1 so that A a = TrNpt^- By using 
even and odd parity combinations of left and right lead 
states, model (1) is reduced to a single-channel Anderson 
model with a resonant level half-width at half-maximum 
given by A = Al + Ar. The spectral function of the 
latter model is required in the calculation of the linear 
conductance, which we describe next. 



III. NRG CALCULATION OF CONDUCTANCE 



II. MODEL 

We consider the simplest model of a correlated quan- 
tum dot, the single-level Anderson model given by the 
Hamiltonian 

H = ^2 £ dv n da ~ 9HbB s d z + Un d ^n d i 

(J 

+ Y^ £k a c\ aa c kaa + X)(*« c L<A + h - c -)- C 1 ) 

kaa kaa 



The linear response electrical conductance G(T, B) of 
(1) is given by'' 1 ' 

G(T,B) = ^ fdu (~j£\ X)T ff («,r,B), (2) 



where 



% (lu, T, B) = 4tt A£ A ^ A a (w, T, B) (3) 

is the transmission function for spin a electrons. It can 
be calculated from the single-particle spectral function 



of the dot A a (w,T,B) = -Jm[G d(T (u) + iS)]/n, where 
Gda{u + iS) = {{d a ;d) a )) is the Fourier transform of the 
retarded single-particle Green function of (1). In Eq. (2), 
e and h are the electronic charge and Planck's constant 
respectively, and /(w) = [1 + exp(/3w)]~ 1 is the Fermi 
function. 

We use the NRG to evaluate the spectral function 
A a {ijj, T, B) via the Lehmann representation 



z 

m.n 

xS(u-(E m -E n )), 



(4) 



where M^ nn are the matrix elements of the spin a lo- 
cal d-electron operator between eigenstates \m) and \n) 
with energies E m and E n and Z — J2 m exp(-/5£ m ) is 
the partition function (see Ref. 34 for details). The usual 
approach is to broaden the discrete spectral function in 
Eq. (4) with Gaussians or Logarithmic Gaussians in or- 
der to obtain a smooth function, 34 ' 41 ' 42 which is then 
substituted into Eq. (2), thereby yielding G(T,B) after 
a numerical integration. A more accurate procedure, in- 
troduced in Ref. 35, is to substitute the discrete repre- 
sentation in Eq. (4) directly into Eq. (2) resulting in the 
expression 



g(t,5) = |eE 



M' 



1 



mn I 



a m,n 



e pE m + e pE n 



(5) 
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r 3 -. This avoids errors from numcr- 



where 7 — ... , , , . 
ical integrations and from an artificial broadening of the 
spectral function and has been shown to give accurate re- 
sults for the conductance. A similar procedure has been 
used in Ref. 43 to extract Ct for the symmetric Ander- 
son model to within 5% accuracy. This uses a full density 
matrix evaluation of the spectral function ■ within the 
complete basis set and is computationally more inten- 
sive than the approach which we use here, whose compu- 
tational complexity is comparable to that of evaluating 
a local static correlation function. 

For the remainder of this paper we shall set the g fac- 
tor g, Bohr magneton /j,b, Planck's constant h, electric 
charge e, and Boltzmanns constant fen to unity, and also 
assume symmetric coupling to the leads (A^ = An = 
A/2,7 — tA). A finite asymmetry A^ ^ A# only influ- 
ences the value of G(0, 0), but not our results for cb and 
ct- Note that in experiments on quantum dots, 1 '' the 
extracted full width at half maximum of the Coulomb 
blockade peaks is given by T = 2 A. 

In evaluating Eq. (5), we used ^-averaging within the 
band discretization scheme of Ref. 36 (see Appendix A) . 
A discretization parameter of A = 4 with n z = 2 val- 
ues for the z-averaging was used and the cut-off for 
the rescaled energies at each NRG iteration was set to 
e c (A = 4) = 30. Figure 2 shows typical examples 
for G(T) versus T/T at B = for a strongly cor- 
related quantum dot (U/A = 16 3> 1) in the Kondo 
(ed = —17/2), mixed valence (Ed = 0) and empty orbital 



Fitting range 


R 2 


c T 


% error 


10" 5 Tb < T < T 


0.818 


0.7447 


87 


lCT 5 To < T < O.lTo 


0.9969 


5.1277 


16 


l(T 5 To < T < 0.05T 


0.99965 


5.8002 


4.7 


l(r 5 T < T < 0.027b 


0.999980 


6.0820 


0.086 


10~ 5 T < T < O.OlTo 


0.9999894 


6.1459 


0.96 


10" 5 T < T < 0.005T 


0.99985 


6.1468 


0.98 


10~ 5 T < T < O.OOlTo 


0.9381 


6.2034 


1.91 



TABLE I. Optimal temperature range for fitting the con- 
ductance G(T,0) to the Fermi liquid form f{T/T ) = a(l - 
c T {T/T ) 2 ) for U/A = 12 and e d = -17/2 using a goodness 
of fit based on the value of R 2 . The latter is defined by 



R 



1 



EIUU-A-T,))" 



^ ?=l{yz - <y>) 2 , where x, = T t /T , Vl = G(T t ,0), 
(y) — 4 5^™ =1 i/i and the number of data points in the fitting 
ranges was n sa 200. The value R 2 = 1 would correspond to a 
perfect fit to the Fermi liquid form. The % error in ct in the 
last column is defined by % error = 100 ■ I— — CT '° xact I where 

J ' C T. exact 

CT.cxact is the exact value at particle-hole symmetry given by 
Eq. (15). From this table we see that the optimal range which 
maximizes R 2 and the accuracy of ct is close to T < 0.02Tb. 
The NRG calculations used A = 4, n z — 2 and an energy 
cut-off e c (A) = 30. 



(ed = (/> A) regimes. The scale To is defined from the 
T = susceptibility of the Anderson model (1) 



X(0) = 1/4T 



(6) 



for all U and Ed- For the case of particle- hole symmetry 
(Ed = —U/2) and strong correlations U ^> 7rA, one also 
has from the Bethe ansatz solution for x(0) an analytic 
expression for T 



T (s d = -U/2) = T (s) a ^/UA/2e 



: ,-TrU/8A+nA/2U 



(7) 



within corrections which are exponentially small in U/ttA 
(see Ref. 13). 

The Kondo scale To is useful in analytic calculations of 
Ct and cb about the Fermi liquid fixed point at T = 0, 
such as those in Sec. IV. With this definition, the mean- 
ing of the coefficients ct and cb is fixed by 



G(T, 5 = 0) _ (T 

G(0,0) ~ CT {Th 

G(T = 0,B) (B_ 

G(0,0) CB \T 



(T « T ), (8) 
(B « To). (9) 



In extracting cb we do not use Eq. (9), but instead use 
the Fermi liquid result in Eq. (17) of Sec. IV. This al- 
lows cb to be obtained directly from the T = occu- 
pancy of the d-level, a quantity that can be calculated 
to high accuracy within the NRG. The coefficient ct 
is extracted numerically by fitting G(T, 0) in the range 
10~ 5 T < T < 2 x 10~ 2 T to the Fermi liquid form 
in Eq. (8) with T as defined in Eq. (6). The range 
T < 0.02Tq was found optimal for this purpose, as we 



now describe. 4 ' Specifically, we fit the NRG results for 
G(T, 0) in the above range to f(x) = a(l — ctx 2 ) where 
x = T/T . We find that a = 2 ± 1CT 5 at the particle- 
hole symmetric point, with a — 2 (in units of e 2 /h) be- 
ing the exact result from the Friedel sum rule. The ef- 
fect of the fitting range on the accuracy of the extracted 
ct and the degree of confidence in the fermi liquid form 
f(x) — a(l — ctx 2 ) may be ascertained quantitatively 
by calculating the R squared coefficient R 2 (also called 
the coefficient of determination and defined in Table I). 
Table I lists R 2 , together with the extracted ct and the 
% error in ct for different fitting ranges. We see that 
the range 10 _5 To < T < 0.02To is close to maximizing 
both R 2 and the accuracy of ct [as compared to the ex- 
act result in Eq. (15) of Sec. IV]. We therefore used this 
range throughout, also for the asymmetric cases. Care 
is needed in the choice of the cut-off e c (A) in order to 
obtain correct results for G(T) in the low-temperature 
limit when using Eq. (5). If e c (A) is chosen to be too 
small, the correct saturation behavior of G(T, 0) in the 
low-temperature limit (i.e. the "leveling-off" of the con- 
ductance) is not obtained. In this case, a fit of G(T, 0) 
to /(T/Tq) shows a drop in R 2 to small values, indi- 
cating a problem. This is remedied by increasing e c (A) 
(the used value e c (A = 4) = 30 was sufficient, whereas 
e c (A = 4) = 12, for example, is not). Thus, the R 2 cri- 
terion can be a useful check on appropriate choices of 
cut-off when evaluating the conductance via Eq. (5). 

In experiments on Kondo correlated quantum dots, To 
is not measurable, and instead one extracts a Kondo 
scale, T£ xp , from the temperature dependence of the 
B = conductance via 



G(T = T° xpt ) 



G(0)/2. 



(10) 



In principle, this T£ xp can be extracted for each gate 
voltage (i.e., for each Ed), but in practice, it is usually 
extracted only at mid-valley (Ed = —U/2) where one is 
sure to be in the Kondo regime for large U/A. This is 
illustrated in Fig. 2 by the dashed lines. 

With this definition of T£ xp , one extracts the experi- 
mentally measured coefficients c™ p and c C g P via 



G[T ' D Q) --l-c^i^LA , (T«T° xpt ), 



G(0,0) 

G(T = 0,B) 
G(0,0) 



Tie 
is 



(ii) 



--1-C 



cxpt 



B 



rriCXpt 



, (B«T° xpt ) 



(12) 



For the particle-hole symmetric Anderson model, the co- 
efficients c™ p and c°g P are related to ct and cb via 



For a precise translation of theoretical calculations of cb 
and ct in terms of To, into experimentally measured 
ones in terms of T£ xpt , one therefore requires the ratio 
rpBxp jj<s a j. m id- valley for all charging energies (U/A), 
which we supply in Sec. VI. 



IV. FERMI LIQUID RESULTS FOR G(T, B) 

For the case of particle-hole symmetry, the coefficient 
ct is known for arbitrary U/A within renormalized per- 
turbation theory about the Fermi liquid fixed point. i,x 
The expression is given by 



ct 



n 4 i + 2(R-iy 



12 



R 2 



(15) 



where R is the Wilson ratio [defined in Eq. (B24)]. In 
the limit of strong correlations, U/A ^$> 1, the Wilson 
ratio approaches 2 and ct takes the well known universal 
Kondo value ct = 7r 4 /16 (see Ref. 49 and 50). In the 
opposite limit U/A — > the Wilson ratio tends to 1 and 
ct acquires the value 7r 4 /12. Evaluation of Eq. (15) for 
general U/A requires knowledge of R, either from Bethe 
ansatz or from NRG. 

Fermi liquid theory allows an exact analytic expression 
for cb to be obtained for all U and Ed- For this purpose we 
use the Friedel sum rule A a (0,B) = sin 2 (nnd<j(B))/TT A, 
where rida(B) is the spin a local level occupancy in a 
small finite magnetic field B <C T at T — 0. Using 
nda{B) — fid/2 + aaB, where rid is the total occupancy 
at B = 0, and the fact that a = | ZMt(g)_pyXg) = % ( ) 
we easily find the exact result (correct to order B 2 ) 
G{0,B) 



G(0,0) 



i-^V(o)s 2 (i 



cot (-g-)) 



CB 



1 - c B 



;(1 



cot ( 



)). 



(16) 



(17) 



16 v v 2 

where x(0) = l/4To has been used. Note that at 
the particle-hole symmetric point (e^ = —U/2), where 
nd = 1, cb takes the universal Kondo value 7r 2 /16 for all 
U. This universal result [obtained also within SPT, see 
Eq. (19)] could be tested in semiconductor quantum dots 
that can be tuned through complete valleys. It would 
then acquire the value 7r 2 /16 at mid- valley for any val- 
ley. The expression in Eq. (17) also shows that cb de- 
creases monotonically with increasing gate voltage away 
from mid-valley, with cb becoming negative on entering 
the mixed valence regime (which we define by the average 
occupation being rid = 0.5). 



cxpt 



c T 



expt 



= C B 



(13) 



(14) 



SPT CALCULATION 



The SPT approach 2, is based on a renormalized per- 
turbation theory on the Keldysh contour 11 '' 2 '"' using dual 



fermions 29 " 31 . This approach can be shown to be cur- 
rent conserving by construction - ' even in the nonlinear 
response regime, as opposed to finite-order perturbation 
theory in the bare parameters J . We compare the re- 
sults of this theory for the linear conductance with NRG 
calculations. The reference system is the interacting 
particle-hole symmetric Anderson model characterized 
by the renormalized Coulomb interaction u = zTq/ttA, 
where z is the wave function renormalization constant, 
and Tq(U) = r^ ; ^(0, 0;0,0) the four-point vertex. In 
order to obtain results for the asymmetric model, an ex- 
pansion in id = z(sd+U/2)/zA = (ed+U/2)/A is carried 
out for the local level Green function up to a given order 



in id and u, currently up to order u i d . As 



with 



increasing U/A such that in the symmetric case u — > 1, 
it follows that u 2 e 2 d increases with growing U/A at fixed 
e.d/A. An outline of the method is presented in appendix 
B, with full details available in the Supplemental Mate- 
rial of Ref.25. The expression for ct, given in Ref. 25, 
and that for cb, derived in Appendix B, are given by 



1 + 2u 2 + (1 - u)(5u - 3)i 2 d 



12 (l + w) 2 (l + (l-u) 2 £~2) 



r,\2sa\2 



1 - 3(1 - u) 



,•^^2^2 



16 (1 + (1 



uy^r 



(18) 
(19) 



where the apostrophe on these indicates that they are 
evaluated by using the susceptibility Kondo scale of 
the reference system (symmetric Anderson model), i.e., 



°T) °B 



c'd are defined via 



(NRG) 

(NRG) 
(SPT) _ 




20 



FIG. 3. (Color online) ct vs U/A for the symmetric Ander- 
son model calculated within NRG (solid lines with symbols) 
and SPT (dashed line). Filled circles show ct using the sus- 



ceptibility scale Tg 



(«) 



while filled squares show ct upon using 
By comparing the value 



the scale from the conductance x£ xpt 

of ct at U/A 3> 1 with the exact one ct = 7r 4 /16 fa 6.i 
we estimate the relative error in the NRG calculation of ct 
to lie below 0.2%, considerably more accurate than previ- 
ous estimates. 43 ''" Inset (a): ratio T£ xpt /T (s) vs U/A. For 
U/A > 1, the ratio T° xpt /T (s) approaches 1.04. NRG param- 
eters were for A = 4 with an energy cut-off e c (A = 4) = 30 
and n, = 2. 



G(T, B = 0) 
G(0,0) 

G(T = 0,B) 
G(0,0) 



=1 - c' 






(T<:T°), (20) 
(B«T *), (21) 



where T = 1/4^(0) is the susceptibility Kondo scale 
for the symmetric model and is given explicitly within 
SPT by Eq. B18 in Sec. B. In order to compare the above 
results with those from NRG, we need to relate the renor- 
malized Coulomb interaction, u, appearing in the former, 
to the bare Coulomb interaction, U , appearing in the 
latter. As outlined in Appendix B, from Eq. (B25) this 
relation is given by u = R — 1 for the symmetric Ander- 
son model, where the Wilson ratio, R, is calculated for 
given U/A from the exact Bethe ansatz expressions for 
the susceptibility and specific heat 53 of the fully inter- 
acting symmetric Anderson model' 1 . Notice that upon 
substituting u = R — 1 into the SPT expression Eq. 18, 
for the particle-hole symmetric limit id = 0, it reduces 
to the exact Fermi liquid result Eq. 15. 




FIG. 4. (Color online) cb (filled circles) and c C g pt (filled 
squares) vs U/A for the symmetric Anderson model calcu- 
lated within NRG NRG parameters were for A = 4 with an 
energy cut-off e c (A = 4) = 30 and n z = 2. SPT (dashed 
curve) also recovers the value cb = 7r 2 /16 for particle-hole 
symmetry. 



VI. RESULTS 



A. Symmetric case (ed 



-U/2) 



Figure 3 shows ct vs U for the symmetric model, both 
in terms of the scale T (ed = —U/2) = Tq and in terms 
of the Kondo scale from the conductance X£ xp (i.e., c^ p , 
discussed below). The former is compared with the cor- 
responding SPT prediction in Eq. (18) and we see very 
good agreement between this and the NRG calculations 
for all U/A. For U/A > 6, the value of ct remains within 
2% of the the universal Kondo value ct — 7r 4 /16 = 6.088. 
The value U/A = 1*3 separates the weakly corre- 
lated (U/ttA < 1) from the strongly correlated regime 
(U/ttA > 1). We see that in the moderately correlated 
regime 6 > U/A > 1, the deviation ct from the Kondo 
value increases, eventually reaching 7% at U/A ps 3. 
For weakly correlated (non-Kondo) quantum dots with 
U/ttA < 1, ct first decreases with decreasing U, reaches 
a minimum at U/A ps 1.7, and then increases to its non- 
interacting value of7r 4 /12 ps 8.117 at U = 0. Note that 
this latter value differs by more than 30% from the Kondo 
value at U ^> A. 

In Fig. 3(a) we show the ratio T£ xpt /T (s) vs U. This 
ratio allows obtaining c^ p from Eq. (13), the coefficient 
measured in experiments, and which we show in Fig. 3. 
Note the very different behavior between c™ p and or for 
charging energies U/A < 5. This is due to the strong de- 
pendence of T£ xp /Tg on U/A in this range of charging 
energies. In particlular c™ p acquires a maximum value 
of « 8.75 at U/A ps 3.5. Since l£ xpt /T W ps 1.041 for 
U/A > 10 [Fig. 3(a)], 4 xpt ps 6.58 for strongly correlated 
quantum dots in the Kondo regime. In contrast, for 
U/A < 10, the scale X£ xp differs appreciably from TJj '. 
Thus, even for nominally Kondo correlated quantum dots 
with U/A ss 4.5, such as those in Ref. I, one finds from 



Fig. 3(a) that X£ xpt /T (s) ps 1.18, so one should expect 
c^ p ps 7.5, which is somewhat larger than the extracted 
value c™ pt ps 5.6 ± 1.2. 

As discussed in Sec. IV, cb is independent of the charg- 
ing energy U for the particle-hole symmetric case, where 
it takes the value 7r 2 /16 ps 0.617, which is also recov- 
ered exactly within SPT [see Eq. (19)]. However, exper- 
iments use the scale X£ xp and measure c^ p as given by 
Eq. (14). This depends on U through the ratio of Kondo 
scales X£ xp /T . For completeness, we therefore show 
the U dependence of c^ pt in Fig 4. For U/A ps 4.5, rel- 
evant for the experiments in Ref. 4, we find c™ p ps 0.89 
significantly smaller than the value c™ p « 5.1 extracted 
from the measurements. As discussed in that paper, the 
large discrepancy between the measured and predicted 
values of c^ p could indicate the importance of the large 
spin-orbit interaction present in the InAs quantum dots 
investigated in Ref. . 




FIG. 5. (Color online) ct vs Ed (in intervals of 0.5A) with 
Ed > —U/2 for several U/A, ranging from strong, [/ > A, 
to weak, [/< A, correlations, and using the scale To. NRG 
parameters were for A = 4 with an energy cut-off e c (A = 4) = 
30 and n z = 2. 




FIG. 6. (Color online) cb vs Ed (in intervals of 0.5A) with 
Ed > —U/2 for several U/A, ranging from strong, [/ > A, 
to weak, [/< A, correlations, and using the scale To. NRG 
parameters were for A = 4 with an energy cut-off e c (A = 4) = 
30 and n z = 2. Inset (a) shows ct/cb vs Ed/ A for selected 
U/A. The dashed line is a guide to the eye and represents 
the universal Kondo value ct/cb = it 2 reached in the limit 
U/A — > oo and particle-hole symmetry. 



B. Asymmetric case (Ed > —U/2) 

1. NRG results 

For completeness, we show the dependence of ct on Ed 
for Ed > —U/2 for U ranging from weakly (U/A <C 1) 
to strongly {U/A ^> 1) correlated quantum dots in Fig- 
ure 5. For strong correlations U/A ^> 1, ct decreases 



monotonically with increasing deviations from the Kondo 
regime, eventually becoming negative after the mixed va- 
lence regime is reached. A similar behavior is seen in the 
local level dependence of cb, shown in Fig. 6. In Fig. 6(a) 
we show the ratio ct/cb versus Sd/A for selected U/A 
which approaches the value it 2 at particle-hole symme- 
try and U/A 3> 1. Notice that for correlated quantum 
dots in the Kondo regime, ct/cb decreases monotoni- 
cally with increasing deviation from particle-hole sym- 
metry. Since Ct/cb is independent of the definition of 
Kondo scale used, it could be a useful quantity to quan- 
tify the degree of correlations in a quantum dot (U/A) 
and the degree of departure from particle-hole symmetry 
for specific gate voltages. 



2. Comparison with SPT 

Figure 7 compares SPT results for the local level de- 
pendence of c' T , as defined in Eqs. (18) and (20), with 
correspondingly defined quantities in NRG. Figure 8 
shows a similar comparison for the quantity c' B defined in 
Eqs. (19) and (21). We see in both cases, that agreement 
between NRG and SPT, holds for i d = [e d + U/2)/A < 
0.25, which is consistent with the SPT calculations car- 
ried out to order u 2 i d . For larger deviations from the 
symmetric point and with increasing Coulomb interac- 
tions, we see an increasing deviation of the SPT results 
from the NRG calculations. In contrast to the NRG cal- 
culation, we also see that the SPT result for c' T ceases to 
decrease monotonically with id for U/A > 3 (correspond- 
ing to a renormalized Coulomb interaction u > 0.76). On 
the other hand, the SPT result for c' B decreases mono- 
tonically with increasing id as in the corresponding NRG 
result. Although we show comparisons also in the region 
ed > 1, by construction the SPT calculation is perturba- 
tive in id and agreement can only be expected in the limit 
id ^C 1, which we find. We also expect that the range of 
agreement between NRG and SPT in both u and id can 
be increased by going to higher-order (see discussion at 
the end of Sec. B4), however, this lies beyond the scope 
of this paper. 



VII. CONCLUSIONS 

In this paper we investigated deviations from the uni- 
versal Kondo scaling in the linear conductance of a corre- 
lated quantum dot due to a finite level asymmetry (i.e., 
deviation of gate voltage from mid-valley) and a finite 
local Coulomb repulsion (i.e., finite charging energy). In 
particular, we determined the behavior of the coefficients 
ct and cb as a function of £d and U within NRG and 
compared these with results from SPT, 25 finding good 
agreement for all U at the symmetric point and reason- 
able agreement for id — {ed + U/2)/A < 0.25 away from 
the symmetric point. Both ct and cb are monotonically 
decreasing functions of the deviation id from the sym- 
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FIG. 7. (Color online) c' T vs Sd/A for several U/A calculated 
within NRG (symbols) and SPT (lines). Legend: column one 
U/A, column two u. c' T is defined in Eq. (18) using the Kondo 
scale in Eq. (20), and the corresponding NRG result uses the 
same Kondo scale for the purposes of this comparison. NRG 
parameters were for A = 4 with an energy cut-off e c (A = 4) = 
30 and n z = 2. 
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FIG. 8. (Color online) cb vs Ed/A for several U/A calculated 
within NRG (symbols) and SPT (lines). Legend: column one 
U/A, column two u. c' B is defined in Eq. (19) using the Kondo 
scale in Eq. (21), and the corresponding NRG result uses the 
same Kondo scale for the purposes of this comparison. NRG 
parameters were for A = 4 with an energy cut-off e c (A = 4) = 
30 and n z = 2. 



metric point id = for all U and an exact Fermi liquid 
expression for cb has been given which is valid for any U 
and Ed- In particular, the coefficients ct and cb become 
negative on entering the mixed valence regime, signaling 
the onset of thermally activated transport which becomes 
pronounced in the empty orbital limit rid ~ 0. 

For the mid-valley conductance, we also determined 



the ratio of the conductance to susceptibility Kondo 
scales T^ xp /Tq , allowing us to relate our results for 

(s) 

ct and Cb in terms of Tq , to the measured coefficients 
c!y? p and c°g P in terms of T£ xp . While for quantum 
dots with U/A > 6, the difference between the two sets 
of coefficients is a constant factor of order unity (e.g., 



cxpt 
C T,B 



Ict.b = (T£ xpt /T (s) ) 2 « 1.08 for U/A » 6), for 
quantum dots with U/A < 6 this difference becomes sig- 
nificant and should be carefully taken into account in 
detailed comparisons of theory with experiment. We ex- 
pect this to be particularly important for semiconducting 
quantum dots since U/A is tunable to smaller values in 
these systems. 
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Appendix A: Alternative derivation of the Campo 
discretization 



In this appendix we give a derivation of the discretiza- 
tion scheme of Ref. 36, following the procedure for general 
energy dependent hybridization functions of Bulla et al. 
in Refs. 37 and 34, which has been used for the NRG 
calculations of the conductance in this paper. 

We start with the single-channel Anderson impurity 
model, given by 

k.a 
k,a 

which may be written in the energy representation as'^ 7 

^ r D + 

H =H mp + V / MeX/X, + a t Mde 

a J-D- 

+ J2 / s( £ ) a L a ^de. (Al) 



Here, a e ^ a and a e , a , obey the standard anticommutation 
relations {a £iCT , a\, a ,} — 5 a . a '5(e — e'), g(e) is the disper- 
sion, h(e) is the hybridization amplitude, and ±D± are 
the upper/lower conduction electron band edges. The 
model (Al) is characterized by the hybridization func- 
tion A(u>) = J^k \Vk\ 2 S(u — £fc)- As shown in Ref. 37, its 



energy dependence may be distributed arbitrarily over 
the functions g(e) and h(e), as long as the following con- 
dition is satisfied 



ah^^M^H) 2 , 

aui 



(A2) 



where e(w) is the inverse function of the dispersion g(e), 
i.e., 

ff(e(w)) = u>. 

Our starting point is the observation by Campo et al. 36 
that a linear discretization of the conduction band with 
a Fourier basis in the discrete intervals leads to a correct 
estimate for A(w), whereas a Fourier decomposition on a 
logarithmic scale as suggested by Krishna-Murthy et al. 
in Ref. 3i systematically underestimates A(w) (or equiva- 
lently the conduction electron density of states p(u) since 
A(ui) = Trp(ui)V 2 for a constant hybridization matrix el- 
ement Vk = V). This underestimation results in an ef- 
fective hybridization function A(u>) — A(u>)/A\ where 
the factor A\ = ^^ i_a-i > 1 is due the discretiza- 
tion and A > 1 is the band discretization parameter. 5 '' 
While this effect may be corrected "manually" for each 
A, it is clearly advantageous to have a built-in procedure 
within the NRG that does this automatically. Campo 
et al. accomplished this within a logarithmic discretiza- 
tion scheme by using a Fourier decomposition in terms 
of non-orthogonal basis functions. As in the case of the 
linear grid, this correctly estimated A(tu). Motivated by 
this, we provide here an alternative derivation of this dis- 
cretization scheme following the procedure of Bulla et al. 
in Refs. 37 and 31 for general A(w). 

We consider the following set of orthonormal Fourier 
functions in each interval of a linear grid, 



i>n, P {ri) = 



e' 2 ^^, if 776 [n,n + l],n = -1,0,1, 



0. 



otherwise. 



The inverse functions are given by ^f np (ri) 
fulfilling the usual orthonormality condition: 



Ipn.piV^n'.p'iv)^] = Sn.n'Sp 



CpW 



(A3) 



We will transform this relation to a logarithmic 
grid such that [n, n + 1] will be transformed to 
D + [A' n - z - 1 ,A- n - z ] for n = 0,1,.... The first in- 
terval [—1,0] is special and transforms to the first log- 
arithmic interval containing the band edge, i.e. to 
D + [A~ Z ,1]. One possible choice, the obvious one, is 
e = D+A-v-z o 77(e) = -\n\e/D+\/\nA - z,n = 
0,1,... (e = D+A- z ^ +1 \n = -1) , but other choices 
are possible for defining the transformation between lin- 
ear and logarithmic grids (and hence 77(e)). Thus, for 



n = 0, 1, . . . we have 



58 



In 



j£j_ 



In 



j± 



I In A 



Vvp(--?^r -*)*«',?'(— j^f-*)* 



In A 



#£ P (.o/c+ 



=c+-*+ ,(e) 



On 7 n' Op,p' 



(A4) 

For the expansion of the negative part of the band we 
use e = —D_A~ V ~ Z and do not reverse the integration 
boundaries yielding the same function inside the integra- 
tion. 

1 ln lf ^TT- 

-VW — 7—r- -«)*n'p'( — r^~ ~«)de 



In A 



In A 



In A 



=<Pn, P (e)/c r , 



=c~-$~ ,(e) 



On z n'Op,p' 



(A5) 

The normalization factor c^ 1 can be distributed freely be- 
tween the new basis <^ p and its inverse 3?„ ]P - a e ,cr is 
expressed in terms of the new basis: 

where 

/+n 
a e . a $+ p (e)de 

K,p,a = J b e ^$~ p (e)de 



where we defined 

D. A-' 



-D_A~ n -' 



Evaluating the anticommutator < a n>Pl<T , a n , , a , > we find 

<\0-n,p,a: a n ' ,p' ,<?' j 

/+n r+n' 

a e ,^t tp (e)de J a\, ^+ / (e')de' 

/+n' r+n 

4,<x' $ n\ P 'V)*/ a e , a $+ p (e)de 

/+n r+n 
J [a e ^al, a ,}^l p (e)^+ p ;(e')dede' 



U7l,7l' ^<7,a' 



1 ^(p-p'X-gj^+z) 



with an analogous expression for < b n ^ P:(J , b' n , , a , >. Set- 
ting {a„,p, (7 ,al ipcr } = {bn.p.a.b'l^^} = 1 fixes the con- 
stants c^ in Eqs. (A4) and (A5), leading to 

|c±| 2 = £ ± A— ^(1-A- 1 ) =dt 
{=D ± {l-K- z ) forn = -l) 



and 



f f I -a A J 1 



In A 
27ri(p-p')+lnA' 



if p = p' 
otherwise, 



with an analogous expression for < b n , P ,cr,b n , , a , >. Thus, 

only for the continuum limit A — > 1 is the above an 
orthonormal basis for all p,p'. M ' However, as we show 
below, an approximate discretized Hamiltonian can be 
formulated in terms of the orthonormal subset of p = 
states only, within which the NRG calculation is carried 
out. We show that for general A(w), (i), only p = states 
couple to the impurity, and, (ii), off-diagonal terms in 
p,p' can always be eliminated from the Hamiltonian by 
a suitable choice of the function 77(e) relating the linear 
to the logarithmic discretization e = ±D±A~ V ^~ Z . 

With the new basis functions we follow the derivation 
of Bulla et al. in Ref. .'-S4, reformulating first the hy- 
bridization part of Eq. (Al) 



'+ 



D. 



+ 11 



h(e)a^ a d€ = ^a n ^^ / h(e)<f)+ p (e)de 

n,p ■* 

r—n 

+ ^2b„ :P ^ / h(e)(f>~ !P (e)de 



The requirement that the hybridization only couples 
to the p = terms can be satisfied by choosing h(e) oc 
*n,o( e ) = 77Jr> which by Eqs. (A4) and (A5) implies 
that p 9^ do not hybridize. Therefore we can choose 
the same h(e) as in Ref. ->4 (i.e. a step function in the 
discrete intervals), 



He) 



,± 2 



1 



±n 



A(uj)du>, 



(A6) 



for D±A~ n ~ z < ±e < D±A~ n ~ 1 ~ z . This choice guaran- 
tees that e(±D±A~ n ~ z ) = ±D±A~ n ~ z [proved by using 
Eqs. (A2) and (A6)] and that the dispersion is linear at 
the grid points g(±D±A- n - z ) = ±D ± A- ,l - z . The first 
part of the hybridization may be written as 



E 



,+£>4 



f 1 V^ 

/ h{e)fla e ,ade = -/=Y\ft (7« a n , , CT + 7« K,o,a) 

J-D- V^ „ 

. /ct/oct 



where 7^ = J A(tu)duj, and the conduction electron 
Wannier orbital at the impurity site is defined as 

/o<7 = J^ J2n 7rt a «,0, CT + ln b n,0,a, 
With Co = En (lit) + (inf ■ 

Next, we reformulate the conduction electron kinetic 



10 



energy term: 



g(e)a\ a a e ^de 



D- 



/+n 



n p,p 



+ b i.p,* b n,p' ,<? I 9( e )i, P ( e )i,p'*( e )* 

/ t / , a n,y,<j a n,p' ,<TSra.p,p' "r" "n,p,er"n,p' ,CTS ra ,p,p' 
n p,p' 






±fi 



5(eXp(e)< p ;( £ )& 



±H 



cjA(cj)- 



rfc( E )i,,,' ( e ) dw 



7r/i(e) 2 
/ ±n A(w)dw 



■-i^i 



I 



±n ^^ e^P-P'X^^+^dLu 

_ (K")|1"A) 2 

J ±n A(co)duj 

For A(cj) = Ao we obtain e(w) = lu and only £ n , = 

ir^V are unequal to zero and agree with the result of 
Campo and Oliveira. Thus, the impurity, which by con- 
struction couples only to the p = state via the hy- 
bridization term, is completely decoupled from the p^O 
states. We now show that the same can be achieved for 
a general A(oj) by a suitable choice of 77(e). 

Following the same derivation as above, but substitut- 
ing 



V(e) 



1,(1) 

rvn 

*4 2) g(z') 



de' 



(A7) 



in Eq. (A3) leads to diagonal £, n ,p,p' for an arbitrary 
A(w). kn and kh, are given by the boundary conditions 
I k W ^y* = n and f tf) ^de = 

n + 1 . c^ and 7^ remain unchanged but 



9(t)<t>i, p (e)<i>t,p>* ( e ) de 



±n J±l.(l) 



d*k y n > 2TTi(p-p') j' k{2) ^de' 



9(e) 

■n ">n U P,P 



de 






£;„ can be obtained by taking the difference of the 
boundary conditions (as defined above), and using 
e(±D ± A- n - z ) = ±D ± A- ,l - z and Eq. (A2), 

, (i) _ I ±n &{u)duj 
" T d±f ±n ^±^ 



As in Ref. 3G the resulting ^ P=P ' — ^ P=P '(C) are given 
by 59 



S n ,p=p'VW ' r±n A(u) 



J ±n A{u)dw 



>-dcu 



(A8) 



The corresponding result, denoted by £, np=p , = 
£« p=p'(B), in the usual logarithmic discretization scheme 
is given by 37,60 



Sn.,p=p' (-") ^ f ± 



J ujA(uj)duj 



(A9) 



J in A(w)dw 

Evaluating (A8-A9) for a flat band with D + = D_ = 1, 
gives 

C,p=p'(C) = ±^(1 + A-^A-'-VAa, n = 0, 1, . . . 

C,p=p'(C) = ±^(1 + (A 2 )- 1 )/^-, „ = -1 

£,„=*>' ( S ) = 4 (1 + A_1 ) A_n ~*' » = 0, 1, . . . 
C,p^(S)=±^(l + (A 2 )- 1 ), n = -l 

We see that £± p=p ,(B)/£± p=p ,(C) is given by the fac- 
tor A\ (A\z for n = —1), indicating that the Campo 
discretization achieves the correct estimate for A(w) via 
a reduction of the effective bandwidth of the discretized 
model. For energies close to the band edge, where the 
Campo discretization gives a different correction to the 
desired one (A^ instead of A&) further corrections are 
needed. <a 

Appendix B: SPT calculation of cb 

We describe the system by the single impurity Ander- 
son Hamiltonian, H = H c + Hd + Hd- C , where we defined 



h c = EE 



ek\C kx<7 C k \a 



\=L,R k,<T 



H d = Y, E d .Jld a + U (d\d t - ^) (d 



(Bl) 



1\ U 



l l a i 



Hd-c= E E [ V k\dlc kXa + V k \cl Xa d a 

\=L,R k,a 



) 



Here, H c is the single-band Hamiltonian for conduction 
electrons at the metallic leads. Hd is the Hamiltonian for 
localized quasi-particle states at the dot, which includes 
Coulomb interaction. Hd- C represents the coupling be- 
tween the dot and the leads. We have defined the spin- 
dependent local energy level Ed a = Ed — crb, with Ed = 
e d + U/2 a small parameter to capture deviations from 
the particle-hole (p-h) symmetric condition td = —U/2, 
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and b = g^sB/2. We build up an SPT calculation start- 
ing from a reference system which is interacting (U 7^ 0), 
particle-hole symmetric (Ed = e^ + [7/2 = 0), and in the 
absence of an external magnetic field (B = 0). 



1. The reference system 

The reference system (Ed = 0, B = 0) self-energy was 
derived in detail in', and is given by the matrix form 



^-la.ui — 






-Y,-+ 



The local Green's function for the reference system is 
given by the matrix 



Scr.u; 



q+~ q++ 



■F(u,T,V)g^ 



(B2) 



with components satisfying 

9--=[l-F(u,T,V)]g r a , 

g-+=-F(u;,T,V)[g r ^-g^J 

gt-=[l-F(0J,T,V)][g^-g a a J 

9t%= -[1 - F(w, T, V)]^ - F(w, T, V)g r a>u . 

Here, the effective local nonequilibrium distribution func- 
tion is shown - ' to be 



F(u,T,V) 



A L f L +A R f R -(i/2)Z- u 



A, 



imsr 



The retarded component of the Green's function is given 
by 

& u = (u> + iA-Tr u )- 1 , 

witl1 a„ u = \SaJ*> and A = A L + A R . It is shown in 
Ref. 25 that the self-energy components for the reference 
system satisfy a similar set of relations as Eq.(B2), in 
particular with 

E-+ = -2ilmYlF(u,T,V) 
E+- = 2amSi[l-F(«,T,V)]. 

The (spin-independent) retarded self-energy of the refer- 
ence system was calculated in detail in , and is given 
by 



K = (i-xa 



->'0 



-iA 



,-2 



X%- 



(- 

\A 



■kT 
~A 



X+- 



C f irTeV Y 



(B3) 



with ££ = (££)*. In Eq.(B3), the "odd" component 
of the spin susceptibility is directly related to the four- 
point vertex T (U) defined in Eq.(B12) by the relation 
X+-(t/)=r /(7rA). 



2. SPT formulation 

A coherent-states path-integral representation of 
the model of Eq.(Bl) in terms of Grassmann fields 
can be obtained following the standard construction. 
On the Keldysh contour, we introduce ipkXa^)^ — 

(pj^r^c+^y and $(t)t = (d-(t),dt(t))\ where the 
indexes ± refer to the time-ordered (-) and anti-time- 
ordered (+) paths, while A = {L,R} labels the two dif- 
ferent leads. 

The resulting non-equilibrium generating functional 
for the model of Eq.(Bl) is 25,62 , 



D[^t ) ^p[4t j ^] e *5[^,^4t,*]_ 



(B4) 



Since the action in Eq.(B4) is Gaussian in the ^^(i), 
^kX<r(t) Grassmann fields, we integrate those in the par- 
tition function Eq.(B4) to obtain, in the frequency-space 
representation, 



V[$i u ,Ku,]e iS & 



(B5) 



In Eq.(B5), we have defined the effective action as 

c+OO 



/J^E*^^*. 



where 



iS u [*L,*™]=iS$ t [*L,* m ,] 

/ + OO J 
-^^|.L(^ + i (A L +A«))a 3 <l CTW (B6) 
-00 Z7r „ 

is the effective action for a particle-hole symmetric (Ed = 
0) and interacting (U ^ 0) system in the absence of an 
external magnetic field (B — 0), and 



iA> 



E 



'k\\ 



U - £feA + iff* 



for A = L, R 



is the coupling with the metallic leads, which in the limit 
of a flat band (p x (u) = p$, 14a = 14) of infinite band- 
width, tends to iT\ — > iirpo\V\\ 2 . In order to construct 
a perturbation theory in the small parameters Ed, B, 
with respect to the reference system defined by the action 
Eq.(B6), let us introduce the dual fermion (Grassmann) 

fields (j>l u = (fo U ,fXi) where, as before, the index =p 
refers to the time-ordered (anti-time-ordered) path along 
the Keldysh contour. We insert the fermionic Hubbard- 
Stratonovich transformation, 29 ' 30 

/| />+00 7 

Vffiau, 4><T(A ex p\ i ^2 7^ [<£L (SauEda^Sa^r 1 

= Dct Ug^Ed^g^r 1 } e-iX.flZ &K„B*,»sK~ 



12 



into the partition function Eq.(B5). Integrating out the 
local fermion held $ CTt j, one hnds that the dual fermion 
bare Green's function is given by 



G /(o) 



-gffw (gffw - E da a 3 ) g au . (B7) 



On the other hand, by functional differentiation of 
the partition function, an exact nonpcrturbative rela- 
tion between the dual fermion dressed Green's func- 
tion GUI = -i 

<JOJ 

E da°3 + {ga.uiEdvdzy 1 G f au) {E da a 3 g a 



au&au) an d the local Green's function 



G'L = -ii&LHu) ^ obtained 



G CT .[ 



(B8) 

The dual fermion Green's function satisfies the matrix 
Dyson equation 



G / = g^°) + G /(0) £ / G f (B9) 

Notice that the zeroth-order solution of Eqs.(B9,B8) is 
Gi°2 = -Edafia + (g^E^&s)- 1 G£°) (E^asg^r 1 
= (g^ci - E da a 3 ) . 
The first-order solution for the Dyson equation (B9) is 

G/« - G / W 4- G f ^T, f G f W 
which upon substitution into Eq.(B8) yields 
GW=Gl°)+G^Gi») 

= [gw- 1 -e/j- 1 + o([e'] s 

= "gW - 1 - S ffiU - £ d(T £ 3 - E£ 

It is clear then that, within this first-order solution of 
the Dyson equation, the perturbed matrix self-energy is 
given by 



S CT , Ba («, B) = (E d + oh) & 3 + S CT]U + E£ . 



3. The dual Fermion self-energy 



(BIO) 



In order to simplify the notation, let us use the multi- 
indexed labels 1 = (wi,o"i,ii) for frequency, spin and 
Keldysh contour %\ = =p indices. Let us define 

D 12 = (2Tr)S UJl - UJ2 5 ai ^ 2 Eda 1 [d'3]i l! i 2 
Ti234= (27r)(5 Wl + W3 _ W2 _ U4 [r o . l0 . a;o . s0 - 4 (wi,W2jW3,a;4)] 



where in this context repeated indices stand for a gen- 
eralized convolution in frequency, spin and Keldysh- 
contour indices. A series expansion of the dual 
fermion self-energy matrix follows from \g — -D -1 ] = 
-D [I - gDY 1 =-D- DgD - DgDgD + . . .: 

E{ 2 = -*r 1234 [gDg] i3 - ^^ 12 34 [gDgDg} 43 + 0(D 3 ) 

(Bll) 

The four-point vertex is given by 

r (o)++++ _ r n A s 

where 

T (U) = C/ + ttA (15 - 3tt 2 /2) (U/nA) 3 + 0(U 5 ). (B12) 

The first term on the right-hand side of Eq.(Bll) pos- 
sesses only two non-vanishing diagonal matrix elements, 

= —iFoEd- <J Z_ a , 
and +iToEd.-o-Z_^ , where we have defined 
'-+ 00 du' 



I duJ ' (\ — 1 



■y—cruj'i •? —ou) ! y—<7D' 



and Z+ + = (Z~~) . Direct calculation of the integral, 
and consistently with the approximation for the reference 
system keeping terms up to 0(Tq) only, we obtain 



E£— (uj,B) = -iT E d> _ a ZZ- 

= -{E d + ab)u\l- 



1 



■kT 
X 



eV 

X 



(B13) 



2l,t2;»3,»4 



Here, V CTia2 . a3(7A {ijj\, w 2 ; W3, W4) is the four-point vertex 
of the reference system. The dual fermion self-energy is 
given by the expression 

E{ 2 = (2Tr)S Ul _ U)2 5 (Tua2 [E£ 1)Wl ] . 



^1234544' [5-D x ] , , g 3 



'3, 



7 /TrTeTA 2 1 

V^J)' 

with the other components given by E£ :++ (a>,B) 
-££-- (w, B), S£+~ (w, B) = E£"+ ( w , B) = 0. 



4. The retarded self-energy 

At the order of approximation of Eq.(BlO) and 
Eq.(B13) , the self-energy components at the local site 
are 

K%M B ) = E ~ - E *> + Zl' ++ ("> B ) 
E--^ (w, B) = E-" ( w , B) + E da + E£~ («, B) 

S+i>,B) = E+-(u;,B) (B14) 

E-+ d ,(«,B) = E-+(w,B). 
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We thus obtain the retarded self-energy from the relation 
£;,*■>,£) = E+i„(w,B) - S++ d (c,B), as follows 

K, Ed (w, B) = (1 - x++)^ + ^d - cr6 - (£ d + erft)« 

eV\ 2 C /TrTeVV 

Xj " 3 ^ A 2 J 

(B15) 



xU- 



-iAte 



,-,2 



X 

2 



ttT 

X 



X 

2 



+ c 



Here, u = zro/(7rA) is the renormalized interaction, 
for z = x++ the wave function renormalization factor 
for the particle-hole symmetric reference system at zero 
magnetic field. The renormalized quasiparticle spectral 
broadening is A = zA. The retarded Green's function 
corresponding to this self-energy is 



= X++{u>-E d + ab + iA-'Z r cru (B) 
Here, we have defined the renormalized self-energy 



tlJB) = -(E d + ab)uh- 1 - 



ttT 

X 



+ 



eV 



A 
+iAu 2 



7 JnTeV\- 



X 



+ c 



eVV C fnTeVV 
Xj ~3 I A 2 J 



In the above, we have systematically included all 
contributions up to second order in the renormalized 
Coulomb interaction u and particle-hole asymmetry id = 
Ed/ A = Ed/ A. This corresponds to approximating the 
dual fermion Green's function by 



G' 



G f 



g£e>g£ 



(B16) 



instead of a selfconsistent solution of the Dyson equation 
G / = G f Q + G^G-f. The self-energy in this equa- 
tion involves a single renormalized four-point vertex, as 
stated by Eq. (Bll). Additional contributions to the self- 
energy are generated by including reducible contributions 
to the four-point vertex. This involves the entire fam- 
ily of "parquet" diagrams. Explicit calculations of these 
higher-order contributions are currently under develop- 
ment, but go beyond the scope of this paper. 



5. Differential conductance G(T, B) 

The differential conductance G(T, B) = dl/dV is ex- 
pressed by the formula 

G(T,B) = ^l'y{- d ^)j:%^T,B) 



Here, the transmission is 

% (u, T, B) = 4n Al A ^ A a (u, T, B) , 
where the spectral function is defined as 

A a (uj,T,B) = --lmGl Ed (m,B). 

IT 

The differential conductance at zero bias, and up to sec- 
ond order in temperature and magnetic field can be cast 
into the form 



G(T,B) , ( T V , (gjXBB/2 "' 

— = 1 - Crp | —^ | - C R | T-V- 






"'^f; (B17) 



Here, the Kondo scale is based on the spin suscepti- 
bility of the particle- hole symmetric system X s (0) — 
( ffW3 ) 2 /4T (s) , with 

X s (0)= ( ^^A d (0)(l + UAd(0) 

2 ttA\ ttAJ 

■£ ttA 
and we have used A d (0) = 1/(ttA). Thus, 



T, 



(») _ 



TTA 



(B18) 



2(1 + u)' 

The coefficients in Eq.(B17)are given by 

, _ 7T 4 1 + 2u 2 + (1 - it) (55 - 3)ij 
Ct ~V2 (1 + h) 2 (1 + (1-w) 2 £ 2 ) 2 

, _ 7T 2 1 - 3(1 - U) 2 £ 2 d 

° B ~ 16(l + (l-u) 2 e 2 )2' 
where we have set Ed = E d /A = (e d + [7/2) /A. 

6. Derivation of the relation u = R — 1 
From Fermi liquid theory, we have the general result 



A d , a (0) = z- 1 A d , a {0) = 



sin (irn da 
ttA 



(B19) 



where z is the wave function renormalization factor, and 
rida is the local level occupancy for spin a. Along with 
this, we have the following Fermi liquid relations for the 
specific heat, spin and charge susceptibilities 

2tt 2 4 7 , ^ 



Xd 



3 

(5Ms ) 2 7 



(B20) 

.l,,„(l»(l -r.1,,,,10)) (B21) 

(B22) 



Xc,d = 2A dta (0) (l - i/I d , CT (0)) 
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where U = z 2 T , A = zA. From Eqs. (B20)-(B22), we 
obtain 



(SMb) 



4 6 

— zXd + Xc,d - ^H^rid, 



n*k% 



(B23) 



and together with the definition of the Wilson ratio, com- 
bined with Eqs. (B20)-(B23), we obtain 



R 



%/"b) 2 Id 



l + EM^O) (B24) 



Substituting Eq.(B19) into Eq.(B24), we obtain 
R=l-\ ~- sin 2 (7rn dCT ). 

7rA 

Let us define u = U/(irA) = zFo/(7rA), with To defined 
by Eq.(B12) as the four-point vertex. Then, we have 



R-l 

sin 2 (jm dcr ) 



(B25) 



Finally, notice that for a particle-hole symmetric system 
n da = 1/2, and hence u = R — 1. 
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